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BOUNDARY-INDUCED DOWNWASH DUE TO LIFT IN A TWO-DIMENSIONAL 

SLOTTED WIND TUNNEL 1 

By S. Katzoff and Raymond L. Barger 


SUMMARY 

A solution has been obtained for the complete 
tunnel-interference flow Jar a lifting cortex in a two- 
dimensional slotted tunnel. Curves are presented for 
the longitudinal distribution oj tunnel-induced down- 
wash angle for various values of the boundary open- 
ness parameter and for various heights of the vortex 
above the tunnel center line. Some (plant dative dis- 
cussion is given oj the use <j these results m calcu- 
lating the tunnel interference for three-dimensional 
wings in rectangular tunnels with closed side walls 
and slotted top and bottom. 

INTRODUCTION 

The problem of determining lift interference 
corrections for slot ted -throat tunnels by using an 
approximate “homogeneous” boundary condition 
has been treated by several authors (refs. 1 to d). 
For both two- and three-dimensional tunnels, re- 
sults have been given for the vertical interference 
velocity at a lifting vortex located on the tunnel 
center line. (See, for example, ref. 1.) 

This report extends the two-dimensional results 
of reference 1 by giving the longitudinal distribu- 
tion of this vertical interference velocity along the 
horizontal line through the vortex, lor various 
vertical positions of the vortex in the tunnel. 
The manifest application of these results is to the 
determination of the tunnel-induced angle of at- 
tack and camber in two-dimensional airfoil testing. 
An important further application is to the deter- 
mination of the tunnel-induced down wash angle 
in three-dimensional-model testing. The feasibil- 
ity of the latter application was pointed out in 


reference 2 and will be further discussed in the 
present report. 

As in the previous tunnel-interference studies, 
the analysis is made for incompressible How. For 
the present problem, which considers only the 
tunnel-induced vertical velocity for given lift, the 
compressibility effect can be taken into a ccou nt 
merely by reading the induced angle at xf\ l — Af 2 
instead of at x (Gothert rule). 

The method of solution of the problem is con- 
sidered to be of particular interest. It uses a gen- 
eralized Fourier series in a form that automatically 
satisfies the boundary conditions. The method 
has apparently not been used in previous wind- 
tunnel interference studies, although similar meth- 
ods of solution have been used, for example, for 
problems involving thermal diffusion (ref. 4) and 
stress diffusion (ref. 5). 

SYMBOLS 

A n} B n constant coefficients in the series expan- 
sion for <p 
l 

,J h 

c chord of two-dimensional airfoil 

c t lift coefficient of two-dimensional air- 

foil 

C area of tunnel cross section 

C L wing lift coefficient 

d slot spacing (distance between center 

lines of slots) 

k height of vortex above tunnel center 

line, expressed as fraction of tunnel 
semiheight 


i Supersedes NAC A Technical Note 4289 by S. Katzotf and Raymond L. Barger, 1958. 
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t unnel seuniheight 

apparent mass associated with flow 
normal to slotted boundary, per unit 

area, of* boundary, ' log;,, esc ^ 
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Mach number 
root of equation (2a) 
root of equation (2b) 
wing area 

tunnel frec-stream velocity 
rectangular coordinates, expressed as 
fraction of tunnel scmihcighl (fig. 1) 
vortex strength, positive clockwise* 
(lift ing) 

t unnel correction factor 
two-diincnsionul-t unnel com»ction fac- 
tor 

that part of tin* total tunnel correction 
factor contributed by the renter row 
of images 

error in tunnel correction factor in- 
curred by using two-dimensional vor- 
tices in place of the* exact images in 
the rows above and below the center 
row; defined only for closed tunnels 
or for tunnels with completely open 
top and bottom boundaries 
tunnel-induced dowmvash angle* 
total perturbation potential (due to 
lifting vortex and tunnel interfer- 
ence) 

potential of lifting vortex 
potential of tunnel interference 
open ratio of slotted boundary (ratio 
of slot width to slot spacing) 


d slot spacing 

<j open ratio of slotted boundary (ratio of 

slot width to slot spacing) 

For the present problem of the* two-dimensional 
tunnel it will be convenient to make all lengths 
nondimensiona! by dividing by the tunnel semi- 
height ft; in other words, tin* tunnel semiheight 
will be taken as unity. A sketch of this tunnel 
showing the coordinate* system and the* lifting 

Tlx* bound- 


(.:-! ) (la.) 

(z -----1) (lb) 

where q— V- 

ft 

T le boundary condition can be* automatically 
satisfied by making use* of a trigonometric Fourier 
series of a more* general form than that commonly 
used. (Mnsiderthe functions sin r u z and cos ll n z 
(//—I. 2, M, . . .) wh(*re the* r„ are the* positive 
roots of 

tan r-\ fjr — 0 (2a) 

and the* ll H are the* positive roots of 

cot //-////-() (2b) 

The ! unct ions sin r tl : and cos R n z form a complete 
ortli igonal set over 1 he interval — 1 ^ zS 1 (n*f. 4). 
The ser»*s 

sin r n z- \ R„ cos R n z 


(bound) vortex is shown in figure* 1 
arv ‘onelition mav be* written 


<p t a $~ ={) 




BASIS OF METHOD OF SOLUTION 

The boundary condition at a slotted wall, 
eonsidered as an e*([iiivalent homogeneous bound- 
ary, is derived in re*fe*re‘nee*s 1. 2, and 3 as 


where 
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1’iuuiK 1. Sketch showing t wo-dimensiomil slolted 
tui i n * 1 of unit semiheight , coordinate axes, and vortex 
lor U ion. 
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satisfies the conditions (la) and (U>) at 3—1 and 
z I, respectively, because 


[s+(j v 1 =X- U rtiu r*+<jr« <* os >V> 

L ' Oz _|,= i n=l 

BJco s A*,,— '/A*,, sin AV) 

[j s ‘ — = X-M~ sin t‘»—U r n <‘ ()S 

+ /Adcos /i n t/A rt sin A w ) 


where the (inantities in parentheses are all zero 
because of conditions (2a) and (2b). Hence, any 
series of the form 


# — XULA,(;r) sin r n z - i B„FJs) cos H a z 

n = 1 


representing a function in the ./*:-plane automati- 
eallv satisfies the boundary conditions (la) and 
(lb). 

It is of interest to note that the arguments of the 
functions of a Fourier-Bessel series can be similarly 


chosen so that each l-erni satisfies t ho boundary 
condition ^+0 ^==0 on the boundary of a circle 


(ref. (V). 


APPLICATION TO THE PROBLEM OF THE LIFTING 
VORTEX 

Boundary conditions. In reference 1 it was 
shown Unit the total perturbation How (vortex 
flow plus boundary interference) in the tunnel 
could be considered to have two parts: (1) an anti- 
symmetric part characterized by a uniform up- 
wash velocity far upstream, an equal uniform 
down wash velocity far downstream, and zero 
down wash along the vertical axis (the 3-axis) 
through the vortex itself, and (2) a flow with uni- 
form downwash velocity throughout the entire 
tunnel, of such magnitude that it cancels the up- 
stream upwash velocity of part (1) and thereby 


X A n sin r n z J rB n cos Ii„z= — 

n ~ 1 


piovides zero net (‘fleet far upstream and uniform 
downwash velocity along tin* ;>axis. 1 he value 
of this uniform downwash velocity (lor the tunnel 

p 

of unit semiheight) was shown to be ^ This 

value was derived in reference 1 for the rase ol the 
vortex at l lie origin; however, the line of reasoning 
is not altered by (‘hanging the vertical location of 
the vortex in the tunnel, so that the same expres- 
sion remains applicable in the general oflcenter 
ease. Accordingly the potential along the 3-axis 

is given bv — .7 tt-Ft J d> ( >ve the vortex loca- 

0 * 4(f/-f- 1 ) 4 

t ion and bv “ .7 , , -below the vortex loca- 

tion. (The potential is here considered to increase 
in the direction of flow, contrary to the usage in 
reference 1; also, a clockwise (lifting) vortex is 
here considered as positive.) 

The problem will be solved for the upstream half 
of t he tunnel: 0, — 1 ^ z ^ 1 ; the solution for 

the downstream half follows from the known sym- 
metry characteristics. For the upstream half, 
then, the problem, as given by the preceding dis- 
cussion, is to find an expression for <p that satisfies 
Laplace’s equation, satisfies equations (la) and 
(lb), approaches zero as .r approaches and 

along the 3 -axis equals l j 4 a ^ )OV(l 

vortex and — . , 1 —7 below the vortex. 

4 (f/+ 1 ) 4 

Solution for total perturbation flow. The total 
perturbation potential is expressed in the form 

<g=X A,/U^sin r H z+Ii H e K * T cos R n z ( 3 ) 

n = 1 

where the v„ and B n are the positive roots of 
equations (2a) and (2b), respectively. The coef- 
ficients A„ and B„ are determined by matching the 
expression with the known potential along the 
~-axis. That is, let 


4(y - 1) + \ 
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Since the terms of the series are orthogonal in the 
range (— 1, 1), each A n is found by multiplying 
both sides of the equation by sin r n z dz and inte- 
grating from — 1 to 1 : 


u J- 1 8in,r -* rfz= “4( fl Ti)J-. 28inr - 3 dz 

r r k . r r 1 

— r n z<h+~\ sin r n zdz 


r f 1 . 

+ 2 I Sin T n Z dz 

Applying the same procedure to the cosine term 
gives 

B n J cos 2 R n z dz = — —J cos R„z dz 

+?],' ™ a "- 2 ,h 




Performing the indicated integrations and simpli- 
fying the results by means of equations (2a) and 
(2b) yields 

a T cos r n k 
" t ~2r n l+g cos 2 r„ 

__ P sin Ii n k 
2R n 1 +(j sin 2 R n 

Substituting these expressions for the coefficients 
in equation (3) completely determines <p in the 
region x<^0. 

Tunnel-interference velocity. The potential 
<p just derived is the sum of the potential = — I- 

2tt 

7 

tun -1 of the lifting vortex and the potential 

<p 2 of the tunnel interference: (p^ (p]+<p 2 . The 

desired tunnel-interference up wash velocity is 

* oz 

found bv subtracting from 

Or oz 


dip 2 dip dip i 

dz dz dz 


■SSjC 


e r,,x cos r n k cos r n z e RnX sin R n k sin h 
1 + <j cos 2 r n 1+^ sin 2 R n 


The entire model may generally be considered 
with adequate accuracy to lie in one horizontal 
plane. Then, since the interference is of interest 
only in this plane, z may be sot equal to k in the 
preceding equation to give 

/ i r " T cos 2 r n k . e R * r sin 2 R n k\ F 
dz 2 jjrf \ 1 + g cos 2 r n "M + g sin 2 J{ n )'2irx 

(x<Q) (4) 

Equation (4) gives the tunnel-induced upwash 
velocity in the upstream part of the tunnel. At 

p 

x— 0, the upwash velocity is — y as already 

mentioned. The tunnel-induced upwash velocity 
in the downstream part of the tunnel follows from 
the symmetry characteristics previously described; 


that is, 



X 

! +(z~ 

k) 2 

(*<<>) 

■8, 

r 

d^ 2 "l 

£>¥>21 

_r_ r i 

4(fif+l) 

'&VJ 

... £>d 

-z.z L 4(J+1)J 

dpi 

1 

r 

c><ftl 

dJ_ 

L~ 

2^+1) 

Sz J-,. f 


when' the symbol J means the value of the deriv- 
ative at the point indicated by the subscript. 

Ca.culated results, — The tunnel-induced verti- 
cal v docities were computed for points along the 
horiz mtal line through the hound vortex, for 
ranges of values of y and k that are considered to 
more than cover the likely ranges of slot designs 
and airfoil locations. The results are plotted in 

figure 2 in nondimensional form as , where 

€ is dowmviish angle and Vis tunnel speed. For 
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Fk;i ! kk 2. Cunt iniieri. 


sional results such as I hose derived in the pre- 
ceding section could be applied in determining 
corrections for three-dimensional models in rec- 
tangular wind tunnels with closed side walls and 
slotted upper and lower boundaries, provided the 
height-width ratio of the tunnel is not too small. 
Specifically, the tunnel-induced angle at any 
point would be the sum of two parts: (1) that dm* 
to the inlinito row of images of the model in the* 
vertical walls, and (2) the two-dimensional cor- 
rections of figure 2, determined as for a two- 
dimensional wing of the same lift at the same 
longitudinal location in the tunnel. The purpose 
of the present section is to review tin* theory of 
this method and to provide some quantitative 
discussion of the method and of the tunnel height- 
width ratios for which it can lx* considered satis- 
factory. For the smaller height-width ratios, for 


which the method is not quite valid, a lurther 
suggestion is given by which the method might 
still be used with an additional modification. 

A complete analysis should include comparisons 
between induced angles determined by this 
suggested method and those determined by some 
exact method similar, say, to those outlined in 
reference 1. Such exact three-dimensional cal- 
culations would he very laborious, however, for 
the case of the tunnel with slotted upper and lower 
boundaries. Accordingly, tin* present approach 
will bo to make exact calculations for the cases of 
completely open and completely closed top and 
bottom boundaries, which are the two limiting 
cases of slotted boundaries, and then to apply 
arguments of plausibility in extending the results 
to the ease of slotted boundaries. 
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Formulation of problem for the two limiting 
cases. -Figure d(a) shows (lie cross section of a 
closed rectangular tunnel (heavy outline) with its 
nearest images. The plus sign in the center of 
the tunnel represents the semi-infinite doublet 
line (or horseshoe vortex of zero span) trailing 
downstream from a point concentration of lift in 
tin* center of the tunnel. The image doublet lines 
are indicated by plus or minus signs according as 
they arc* the same as or the* reverse of the center 
one. The tunnel-induced downwash or upwash 
velocity at any point in the tunnel (sav at the 
lilting element itself or at some downstream point 
representing the tail location) is the net induced 
velocity due to all the images. The concepts 
underlying this figure and the application of this 
theory have been widely used in wind-tunnel- 
correction analysis (for example, refs. 7 and X) and 
will not be further discussed here. 

As pointed out in references 7 and S, tin* field 
of a doubled line* at a point sufficiently distant 
from it is practically the same as the field of a 
narrow horseshoe vortex having the same moment 
(same product of vortex strength and span- 
in other words, the same* lift). Suppose, then, it 
were 4 valid to replace all the image doublet linens 
in the' rows above and below the center row by 
horseshoe vortices of span equal to the* tunnel 
width. The trailing vortices would all cancel in 
pairs and only the bound vortices, now joined to 
form continuous infinitely long vortices, would 
remain, as indicated in figure d(b). This set of 
joint'd hound vortices, however, will be recognized 
as identical with the' set of image hound vortices 
that arc* used to calculate corrections for a closed 
two-dimensional tunnel (ref. 9), If the* image 
system of figure o(b) were indeed an accurate 
approximation for the image* system of figure* 
o(a), it would follow that the tunnel-induced 
angle* for this ease is the sum of (1) the effect of 
the* image's in the center row and (2) the two- 
dimensiona] correction, just as suggested in the 
opening paragraph of the preceding section. 

The corresponding image diagrams for the tun- 
nel with open top and bottom are shown in figures 
4(a) and 4(b). The argument is essentially the 
same as for the closed tunnel and will not be re- 
peated. 

The problem now is to determine* for what 
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(a) Exact image system. 

(b) A])]) -miniate image system. Doublet, lines above* and 
below the* center have been replaced by horseshoe vor- 
tices < f span equal to the tunnel width. 

I'hu hk 1. — Image systems for a doublet line in a closed 
rectangular tunnel. 

height- vidth ratios the image* systems of figures 
d(b) aid 4(b) satisfactorily approximate the 
exact image systems of figure's ;t(a) and 4(a), 
respect voly, so that, the suggested enlemlation 
procedi.re would he accurate. 

Calculations of the error. In order to investi- 
gate th s problem, calculations were* made of the 
error in the induced-angle factor 8 along the tun- 
nel con.er line that would result from using the* 
image systems of figure's ;t(b) and 4(b) instead of 
those of figure's d(a) and 4(a), Respectively. 
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(b) Approximate image system. Doublet lines above and 
below the center have been replaced by horseshoe vor- 
tices of span equal to the tunnel width. 


Ficurk 4 . —Image systems for a doublet, line in a rec- 
tangular tunnel with closed side walls and open top 
and bottom. 

Thu calculations wore made for tunnel height- 
widtli ratios of ().(>, 0.75, 1.0, and 1.2. These* 
calculations for the zero-span horseshoe vertices 
(or point concentrations of lift in the center of 
the tunnel), however, represent the least favor- 
able case. In order to investigate the problem 
for conditions more nearly representative of typi- 
cal wing tests, calculations were also made for 
horseshoe vortices of span equal to 0.4 the tunnel 
width. As a somewhat extreme case, some 


calculations were also made for horseshoe vor- 
tices of span equal to 0.8 the tunnel width. 
For tin* limiting case, for winch the horseshoe 
vortex lias a span equal to the tunnel width, the 
exact image system would he identical with the 
two-dimensional image system, and thus the 
error would be zero for any height-width ratio. 

The calculated errors arc plotted in figure 5 
as the error A5 in the usual tunnel correction 
factor <5 in the equation 

.s r 

l 

Specifically, A<5 is the* induced down wash factor 
given by the exact image system minus the in- 
duced downwash factor given by tin* approximate 
image system. The approximate image system 
gives too small a downwash angle (too small a 
factor) for the open top and bottom and too large* 
a downwash angle (too large a factor) for the 
closed tunnel. In other words, A8 is positive for 
the tunnel with open top and bottom, and is 
negative for the closed tunnel. The value of 
A5 at *=() represents the error at the lifting 
element itself; values at positive values of x 
represent the error in downwash behind the wing. 
The upstream (negative x) parts of the curves 
are not shown, since the curves all have the same 
type of symmetry as the curves of figure 2. 

As an indication of the order of magnitude of 
the error that might arise, consider the largest 
value of A<5 in figure 5, 0.092, for the tunnel with 
a height-width ratio of 0.0. If tin* wing area 
8 is one-tenth the tunnel cross-sectional area ( r 
and the lift coefficient ( \ is 1.0, the corresponding 
error in tunnel-induced angle is 0.092 X 0.1 XI. OX 
57.;*° — 0.52°. The error decreases rapidly with 
increasing height -width ratio and with increasing 
span of the horseshoe vortex, and it is well within 
usual experimental accuracy for height -width 
ratios of about 1.0. 

Application to the slotted tunnel.- -For those 
tunnel height-width ratios for which the results 
show no practically significant differences be- 
tween the field of the true image system and the 
field of the system of infinite line vortices that 
approximates it, it would seem, reasonable to 
suppose that the effect of slotted top and bottom 
boundaries would be similarly two-dimensional. 
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In those eases in which the difference does become 
practically significant (for example, for tunnel 
height-width ratios less than 1.0 with the larger 
lift coefficients and wing areas) the matter is 
less obvious. A suggested approach is as follows: 

Let the curves of figure 5 be considered as 
correction curves rut her than as error curves. 
That is, suppose that the tunnel-induced angle 
(for the completely open or completely closed 
top and bottom) is considered to be the sum of 
the following three components: 

(1) That due to the images in the center row 

(2) The two-dimensional induced angles cor- 
responding to the infinite line vortices of figure 
‘1(b) or 4(b) 

(3) The correction indicated by figure 5, 
which, essentially, corrects the two-dimensional 
values of component (2) to the correct three- 
dimensional value corresponding to figure 3(a) 
or 4(a) . 

Then a parallel procedure for the tunnel with 
slotted top and bottom would bo to suppose the 
tunnel-induced angle to be the sum of the follow- 
ing three components : 

(1) That due to the images in the center row 

(2) The two-dimensional induced angles given 
by figure 2 for the appropriate value of (j 

(3) A correction determined by interpolating 
between the closed and open eases of figure 5. 
The interpolation could be according to the 
“degree of openness” of the slotted boundary as 

indicated, saw bv the factor that deter- 

£/+l 

mines the two-dimensional correction at jr=0. 
This factor is zero for the (dosed tunnel (g =co ) 
and 1 for the open tunnel (r/— 0). According to 
this suggested interpolation procedure, the con- 
tribution of this third component, to <5 is 


or 

If the slot design is such that <j— 1, this third 
component will be approximately zero. 

In order to help in the application of this pro- 
cedure, the first component 8 n due to the images 
in tin* center row, has been plotted in figure (>. 
The net induced down wash factor, finally, is 

5:2 i "i ' ^&open 
1 T <f 

where 8 C is given by figure (>, is given by 

figure 2, and A 8 0V(n is given by figure f> (positive 
values). 

Asymmetric locations of the model in the tun- 
nel. A sting-mounted model might, at high 
angles of attack, be considerably above the tunnel 
(•(Miter line. Figure 7(a) shows the image system 
for such an asymmetric case, for a point concen- 
tration of lift in a closed tunnel. The correspond- 
ing approximation by means of infinite vortices is 
shown in figure 7(b). No calculations wore made 
for tins ease, but it is apparent from the image 
system that the error is very nearly the average* 
of that for a tunnel whose height is the distance 
m\ to the nearer row and that for a tunnel whose 
height is the distance m 2 to the farther row. 
Except for this one modification, the discussion 
for the slotted-boundary tunnel follows as before. 

Yawing a sting-mounted model also tends to 
increase the error. Figure 8(a) shows the image 
system for an extreme case in which the point 
concentration of lift has been brought close to one 
of the side walls. The image system in this case 
is very nearly the same as for a centered model in 
a tunnel of half the height-width ratio (fig. 8(b)), 


A8 r 


■A( 


A8 oprn — A8 C 


d ) 


where A 8 c i osed and A8 opcn are the values of A5 given 
by figure 5 for the closed and open cases, respec- 
tively. Furthermore, since, as figure /> shows, 
A 8 open and A8 c i 0U(f art* almost equal, except for 
sign, this expression may be approximated as 


-A5, 


opcn~\~ j open 



Longitudinal distance from lifting element, tunnel widths 


Frei io; 6. — Contribution 8 C of the center row of doublet 
images lo the tunnel correction factor. 
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(a) 


+ 




(a) Doublel near the* side wall. 

(b) Don!) <‘t at center of a tunnel having twice the width 
of the tunnel diagrammed above. 

Figure 8 -Image system for a doublet near the sidewall 
of a cl< sod rectangular tunnel, and comparable image 
system for a doublet at the center of a tunnel twice as 
wide as t he first . 


(b) 


(a) Exact image system. 

(b) Approximate image system. Doublet lines above and 
below the center have 1 been replaced by horseshoe 
vortices of span equal to the tunnel width. 

Figure 7. Image systems tor a doublet line above* the 
center of a closed rectangular tunnel. 

for which the error is grout ly increased. Fortu- 
nately, lateral displacement of the lifting point 
to the extent shown in figure 8(a) is far in excess 
of what might he considered a realistic displace- 
ment due to yawing a model in a wind tunnel. 


CONCLUDING REMARKS 

A practical mathematical method is described 
for dote ■mining tin* flow in a two-dimensional 
slotted wind tunnel containing a lifting wing. 
Calculat 'd results are given for the longitudinal 
distribution of tunnel-induced angles for various 
slot geometries and for various vertical locations 
of the ving in the tunnel. Some 1 quantitative 
discussicn is given of the application of these 
results t ) the testing of three-dimensional lifting 
models n a tunnel with closed side walls and 
slotted top and bottom. 

Lang lev Research Center, 
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